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Abstract 

We show how reaction coordinate path lengths affect the relax- 
ation efficiency of a complex system. To this purpose, we consider the 
metric contributions to the transition rates. These metric contribu- 
tions preserve informations about the geometry change of the system 
at the barrier crossing and, therefore, are directly related to the path 
lengths. The output of the present work can lead to identify a way 
to enrich our knowledge about the ability of the complex systems to 
relax only to preferential structures. 

The different aspects of the topographies and topologies of multidimen- 
sional potential surfaces of polyatomic molecules, clusters, and nanoscale 
particles and the way they govern the dynamics and phaselike behavior of 
these systems are currently receiving a great interest; for reviews, see |^, § 
and the references therein. As a result, a new picture of their energy land- 
scape is emerging which greatly enhances our understanding in this pressing 
issue. Despite the enormous power to do computations and many methods 
for extracting useful information, the problem of relating potential surfaces 
of complex systems to their relaxatio is far from a complete solution. 



1 



For complex systems, development of better potential functions still re- 
mains clearly a priority. We can say that, at least, starting from available ex- 
pressions for the interaction between and among the system constituents, ac- 
curate potential energy surfaces (PES) may be constructed for many systems. 
Finding minima and saddles on potential surfaces has become a straightfor- 
ward matter by employing some of the most widely used algorithms or com- 
binations thereof. 0, ^ |^ As it evidenced the problem shifted from how 
to find minima efficiently to what method should be employed to provide a 
better connection between the topography and the dynamics on the surface. 
Statistical methods, [|, ^ |T0[ especially that based on the topography diagno- 
sis procedure, [H achieved considerable success, caying the study of complex 
systems to the point where one can begin to infer dynamics of flow on the 
surface from the general characteristics of the topography. A serious short- 
coming of the existing dynamical studies, however, has to be pointed out: 
virtually, all of them describe the thermally activated barrier crossing by the 
traditional transition state theory (TST) rate model which deals only with 
the height of the barrier and the densities of initial and saddle states. No 
information about the geometrical structures of the initial and saddle sites 
in this approach. From this point of view, there is a limit to the precision 
with which we can characterize the capability of complex systems to relax 
preferentially to one of only a limited number of geometrical structures from 
the vastly larger variety that the system might exhibit. However, there are 
other ways to compensate for this lack of information. It was shown, for 
example, that the topographical properties of the potential energy surface 
(PES) determine the extent to which the system is either glass-forming or 



focusing on regular structures. P, K 12] In addition, qualitative interpreta- 



tions of the interplay between structure and dynamics have been made of 
the path length distributions between connected minima along the reaction 



coordinate. |T^, |T3| 



For a better understanding of dynamics, it is now useful to investigate 
the relation between the variety of path lengths between neighboring minima 
along the reaction coordinate of a specific PES and the fiow of probability 
between those minima. We attempt in the following to get more insight into 
this particular aspect by invoking a generalization of transition-state theory 
(TST) to express the transition rates on the PES.[|I5 To this purpose, we 
include the metric contributions to the transition rates and speculate on their 
relation to the reaction coordinate paths brought about in this way. We work 



2 



out a simple three- level system with a model surface based on a "nearest- 
neighbor" connection pattern and focus on the ability of the system to find 
the global minimum under the assumption of unequal path lengths between 
connected minima. 

Let X denote the reaction coordinate of a complex system with N degrees 
of freedom whose dynamics is governed by an appropriate potential surface. 
The reaction coordinate x is a function of all the degrees of freedom x = 
(gi, ...qN',Pi, ■■■Pn), where qi stands for the i*^ configuration coordinate and 
Pi for its momentum. In the following, we restrict our considerations to a 
small region of the PES containing two states of local stability Oi and O2, 
for which the reaction coordinate takes the particular values Xi and X2 with 
the corresponding energy levels Ei and E2, respectively. In general, these 
two domains of attraction might be separated by barriers containing one 
or possibly more saddle points, possible unstable limit cycles, or even more 
complex unstable attractors, including combinations thereof. For most of 
the following we shall restrict our discussion of the theory to the situation in 
which two adjacent attracting basins are separated by a single saddle point. 
In the present case, we assume that between the local minima 61 and 62 
there is a transition state at X(i_2) = Xa- 

Generally, the forward TST rate wi-2 is given byp!5[ 



{6{x - Xa) X {a)e[x (a)]) 

"^-^ {e{xa-x)) ' 

where 6 and 6 are the usual delta- and step-functions. Here, the average (...) 
denotes an equilibrium average over the canonical probability density. The 
integration over the momenta is straightforward and leads to 

/o^/^^-l/2 jx-Xg) |Vq (x)|] 

= ^'"^^ [e{xa-x)] ' 

where |Vq(x)|^ = Y.i(-0r) Qi ^'^^ mass-weighted coordinates 



{Qi = <lii^y^^- [•••] indicates an average over the coordinates only. Fur- 
ther simplifications in (2) can be achieved by integrating over the coordi- 
nates of the center-of-mass position R and all orientations Q relative to a 
reference configuration of the system. This follows the coordinate transform- 
ing {Qi, ...jQn) {R,^i fi---fm), where f are the internal coordinates (in 
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number, m = N — 6). The internal coordinates allow us to describe forces 
acting in the system in terms of the potential function U (f). For example, 
the number of particle in the well [9 {xa — x)] can be written after suitable 
integrations as 

[e{xa-x)] = Jdh...dUJ{i)9{x{i))exp(^-^^ . (3) 

J(f) is the integrated Jacobian of the coordinate transformation and com- 
prises the determinant of the metric tensor times the volume of the system 
and factors resulting from integrations over the orientational degrees of free- 
dom. 

In the low temperature limit {PU (/) ^ 1, with j3 = ■^^), eq. (2) can 
be evaluated by a Gaussian steepest-descent approximation. The procedure 
needs to expand the effective potential {7e//(f) = [/ (f) — kBT\nJ{i) to 
second order at the saddle point (a) and at the local minimum (6i) . After 
some simple algebra, the final formula for the TST forward rate can be given 
in the form 



Wl-2 



^^^^Vvwexp(-/5i?i_„) . (4) 



27r VJi/ni^l'A 



Note that in the case of multiple transitions states we have to count each 
contribution separately and write the rate as a sum over pathways. H The ~ 
sign is used instead of the equality sign because of terms of order O Q^^ com- 
pared to unity have been neglected in above. M is the mass matrix and the 
indices (a) and (1) indicate that the corresponding quantities are evaluated 
at the saddle point and the local minimum 1, respectively. Ei_a measures 
the barrier height for the forward transition. A, are the eigen values of the 
force constant matrix, that is the covariant second derivative of the energy. 
The Jacobians Ji^a in eq. (4) contain information about the volume, bond 
lengths and orientations of the system at the local minima Gi and the sad- 
dle point, respectively. Equivalently said, the metric contributions Ji^a are 
directly related to the values of the reaction coordinate at these particular 
sites the system is visiting in the evolution through the configuration space. 
Under these circumstances, the above contributions modify the rate constant 
by a multiplicative factor. This factor is the ratio between the two Jacobians 



4 



corresponding to the atomic arrangements at the saddle point and at the po- 
sition of the initial well, respectively, all taken in the same reference frame. 
Within the Gaussian steepest-descent approximation for integrations in the 
phase space, this Jacobian ratio can be taken constant, for any set of coordi- 
nates employed. The backward rate may be obtained straightforwardly 
by replacing J2 and i?2-a for Ji and E'l-a; respectively. 

In the following we shall apply the TST rate formula to a simple landscape 
case in which the two states of local stability Gi and ©2 are connected to the 
global minimum, 63. The global minimum has the energy [E3 < -E'1,2) and 
is placed at along the reaction coordinate. A barrier placed at X(^2-3) = 
Xb separates the global minimum from the nearest-neighbor attractor ©2. 
All the considerations in above apply equally to the forward {W2-3) and 
backward (103-2) transitions between the attractors ©2 and ©3. These can 
be constructed in a similar manner by using the parameters Jj,, J2,3, '^'^^ 
and £'2-6,3-6- 

The general procedure describing the dynamics on the PES relies in prac- 
tice on the explicit knowledge of a master equation governing the time de- 
pendence for the single-event probability Pj (xj, t) of the reaction coordinate. 
For the PES described above, the system of master equations has the form 

dPi 



dt 
dt 



= W2-1P2 {X2) - W1-2P1 {xi) 

= - {w2-\ + W2-3) P2 {X2) + W1-2P1 (xi) W3-2P3 (a^a) (5) 

^ = W2-3P2 {X2) - W3-2P3 {X3) . 

The practical problem is to determine the extent to which the dynamics of 

the system is affected by the metric contributions J to the corresponding 
TST rates Wi^j = 1, 3^ . To do so , we have to solve the above system 
of kinetic equations under certain initial conditions. For simplicity, let us 
assume that at the initial moment {t — 0) the system is in the domain of 
attraction ©1 with the probabihty Pi {t — 0) — 1. Here, the evolution of the 
reaction coordinate through the configuration space starts with the highest 
probability from x = Xi. 

For the present case, analytic solutions of eqs. (5) are available: 



Pi{t) = a 



1+ *£L^f^V«._^e-- 
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1 + 



-012* 



X 



(6) 



1 + 



252 



(2 + x) e-"i*-(2 + ^) e- 



which were obtained in terms of eigenvalues and eigenvectors of the charac- 
teristic system of equations. The (nontrivial) eigenvalues are given by 011^2 = 

^ {Si ± ^2) with -Si = 1 + ^1 +52 + gz, S2 = ^/ (1 + 91-92 -93^ + 4.93, 

_ wi-2 _ W3-2 _ W2^_ p^Yl the other constants entering eq. 

W2-1 ' W2-1 ^ 102-1 " ^ 

(7) are given by x = 5*4 - 5*2, = S'2 + 5*3, ip ^ S2 - S3, tp ^ S4 + S2 and 



a 
S5 



So 



9i 



S^ — So and 5*6 



S:i—S2 \ S4—S3 
S,i 
52 



91 



C2 



S3—S2 

S4—S2 , 



with S\ 



3,4 



1 T 51 ± ^2 + ^3, 



The energy levels corresponding to the minima ©i, ©2 and ©3 of the stan- 
dard potential energy surface employed in the present study are disposed as 
in a steeper funnel obeying the following sequence Ei > E2 > E3. Each 
transition state connecting two adjacent minima lies an energy 1^ above the 
nearest uphill minimum. Therefore, one sets £'2-0 ~ £"1-0 = h > and 
£'3_fe — i?2-a = ^2 > 0, respectively. In addition, we assume in this model 
that the potential energy barrier opposing the escape from the global mini- 
mum is higher than that opposing the escape from the minimum ©2 towards 
the minimum ©i, (£^3-6 > E2-a)- It is also assumed that the relaxation rate 
from ©2 towards the global minimum is energetically greater than the up- 
hill escape -£"2-0) (-E'2-a — -£"2-6 = ^3 > 0). In our computation A*^^^'^) denote 
the vibrational frequencies of the reaction coordinate in the corresponding 
minima and we shall keep them at constant values. As for the mean vibra- 
tional frequency of the transition state, we make the usual assumption, [?] 
that is the geometric mean of the vibrational frequencies of the two minima 
it connects, A('^) = VA(i)A(2) and A^^) = v^WW), respectively 

By using probability distributions Pi {i = 1,2,3) as derived above, we 
have investigated the tendency of relaxation of the complex system under 
various circumstances for reaction coordinate paths. In simulating the dis- 
tribution of path lengths we speculate on the direct relation between the 
Jacobian J (f ) and the reaction coordinate, as can be seen below. The re- 
sults are displayed in Fig.l. 

Going back to the purpose of the present paper, we focus on the role the 
metric contribution plays in the relaxation behavior of the system and assign 
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to the Jacobians Ji, J2, J3, J a and Jh the following sequence of numbers 
Ji : J2 : </3 = 1 : 2 : 3, Ja ■ Jb = 4: : 5. Roughly, these numbers should cor- 
respond to a distribution of path lengths along the reaction coordinate with 
increasing step sizes towards the global minimum. The pictorial correspon- 
dence in nuclear configuration space is therefore characterized by sizeable 
rearrangements of the system components on their way to relax towards the 
ground state. We are interested in monitoring the relaxation dynamics of 
the system on the particular energy landscape as described in above. The 
values of the parameters /i_3 have been tuned to achieve a rapid satura- 
tion effect on a scale of 100. The time evolution of the probability P3 for 
which the system relaxes to the global minimum is displayed in Fig. 1 (see 
the curve labeled (1)). One can observe that the increase in the popula- 
tion of the ground state is sharp and the probability reaches a plateau. We 
now modify the metric contributions to the relaxation rates by considering 
Ji = J2 = J3 = AJa = BJf, = 1 which corresponds to an equally spaced dis- 
tribution of minima along the reaction coordinate. All the other parameters 
in above remain at constant values. Looking at Fig. 1 (see the curve 2), we 
observe that the accumulation in the global minimum is much slower in this 
case. The time evolution of the probability P3 does not reach saturation on 
the same time scale. This indicates that the backward rates of escape from 
attraction domains are higher for the present case in comparison with the 
former and, the system spends considerably more time now moving uphill 
on the PES. The geometrical structures corresponding to minima Gi and G2 
attract the system at rates comparable to that of the ground state structure 
even if the latter is energetically more favorable. Note that, the balance may 
decisively be turned around if, by chance, the numbers of pathways towards 
61 and 02, respectively, are large enough to compensate the energy gaps 
by entropic contributions.]^ The efficiency of relaxation towards the global 
minimum can be reduced dramatically by inverting the numbers in the se- 
quence Ji : J2 '■ J3 from those assumed in the first example. This becomes 
Ji : J2 : J3 = 3 : 2 : 1, and the pictorial correspondence of the PES may be 
that of a steeper funnel with a wide step at the top, between the minima Bi 
and 02 and a narrow one at the bottom, between the minima G2 and 63. 
(We assume : = I, for the present case). As can be seen in Fig.l (see 
the curve 3), the probability P3 of accumulation in the global minimum is 
even worse than that of the previous case. In turn, a narrow barrier between 
the global minimum 63 and adjacent local minimum G2 combined with a 
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broader one between B2 and G3 should result in a longer survival of the 
system in the well around energy level E2. This is demonstrated by curve 
3' in Fig.l which shows that the accumulation in the domain of attraction 
02 is still rising over the entire time scale for the above values of the metric 
contributions. This behavior contrasts with the two previous situations (see 
curves 1' and 2' in Fig.l) where, after an initial increase in population, the 
attractor 62 starts, more or less suddenly, to depopulate. 

In conclusion, we can say that the efficiency of relaxation towards the 
global minimum can be much affected by the metric contributions to the 
transition rates. The metric contributions are directly related to the dis- 
tribution of the path lengths by the integrated Jacobians. Therefore, these 
play, in concert with the barrier heights, an implicit role in classifying the 
archetypal energy landscapes. [0, 0, ||, |12[ In addition, we can say that in- 



corporating the metric contribution to the TST rate is the appropriate way 
to combine the height of the barrier and the densities of initial and saddle 
states with the differential path length corresponding to the transition of 
the system between these specific stationary points. The present approach 
allows the TST rate to be implicitly related to the geometrical structures in- 
volved in the min-saddle-min transition of the system. We intend to use the 
information so obtained to identify a way to enrich the precision with which 
we can characterize the ability of complex systems to relax, preferentially, to 
only a limited number of geometrical structures from the vastly larger variety 
that the system might exhibit. 
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Figure captions 

Fig. 1 - Time evolution of the probability P3 (curves 1, 2 and 3) and 
P2 (curves 1', 2' and 3') for three different sequences of path lengths be- 
tween connected minima. The initial population in the highest minimum 
was assumed equal to unity. Pi (t = 0) = 1. For explanation, see the text. 
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